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Constructions and Noise Threshold of Hyperbolic 

Surface Codes 

Nikolas P. Breuckmann and Barbara M. Terhal 


Abstract —We show how to obtain concrete constructions of 
homological quantum codes based on tilings of 2D surfaces 
with constant negative curvature (hyperbolic surfaces). This 
construction results in two-dimensional quantum codes whose 
tradeoff of encoding rate versus protection is more favorable 
than for the surface code. These surface codes would require 
variable length connections between qubits, as determined by 
the hyperbolic geometry. We provide numerical estimates of the 
value of the noise threshold and logical error probability of these 
codes against independent X or Z noise, assuming noise-free 
error correction. 

Keywords—Quantum error correction, hyperbolic surfaces, sur¬ 
face code 


I. Introduction 

An essential component of scalable quantum computing 
is the ability to protect quantum information reliably from 
decoherence. Such protection can in principle be achieved 
via the use of quantum error correcting codes, assuming a 
sufficiently low physical noise rate and a sufficiently large 
encoding overhead HI. A particularly interesting coding ar¬ 
chitecture is that of the surface code or the toric code first 
formulated by Kitaev |j2l, d, d, Q. The toric code is an 
example of a quantum low-density parity check (LDPC) code, 
meaning that the parity check operators used in quantum error 
correction only involve a constant number of qubits and each 
qubit participates only in a constant number of parity check 
measurements. 

The ideas introduced by Kitaev have been generalized to a 
construction in which a tiled, closed, n-dimensional manifold 
M is translated into a quantum LDPC stabilizer code such that 
all parity checks are local on M. Codes of this type are called 
homological stabilizer codes (see Sectionlllli. The surface code, 
based on tiling a euclidean surface, is an attractive code to 
realize experimentally due to its high noise threshold of almost 
1% and the planar qubit layout and inter-qubit coupling (see 
e.g. recent work in superconducting qubits aimed at making 
elements of a surface code architecture ||6l, jT], IfSlI). 

However, the surface code also has several non-optimal 
features which directly result from the fact that it is based 
on tiling a 2D euclidean manifold. In ||9l it was shown that 
any 2D stabilizer code with n physical qubits has distance 
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Fig. 1; The {5,4}-tiled hyperbolic plane in the Poincare 
disc model. Each tile is a pentagon and four tiles meet at a 
vertex. All lines meet at an angle of tt/2. 


d < 0{y/n) and in ifTOl it was shown that for any [[n, 

2D stabilizer code the number of encoded (logical) qubits k 
and the number of physical qubits n and the distance d obey 
the trade-off relation kd? < cn for some constant c. Both these 
results assume that physical qubits are laid out according to an 
euclidean geometry. Note that these sets of results do not only 
pertain to homological stabilizer codes but hold for any LDPC 
code with a planar qubit lay-out with 0(1) distance between 
qubits. 

This rate-distance trade-off bound is achieved by the surface 
code in which one can, for example, encode k qubits into k 
separate surface code sheets, each with d^ qubits and distance 
d, leading to a total number of n = kd^ qubits. Clearly, the 
encoding rate k/n approaches zero when one tries to encode 
better qubits with growing distance d. 

In ini Fetaya showed that two-dimensional homological 
codes, based on tilings of arbitrary two-dimensional surfaces, 
obey the square-root bound on the distance, i.e. d < 0{^Jn). 
In m it was shown that with a family of homological codes 
based on 4-dimensional manifolds with non-zero curvature, 
one can go beyond this square-root bound and encode a qubit 
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with distance scaling as ^ (nlogn)^/^. In addition, ifT^ 
showed that there exist hyperbolic surface codes which have a 
constant rate k/n> ci while each encoded qubit has distance 
scaling as C 2 logn for some constants ci and C 2 . Delfosse ifTSl 
has shown that the logarithmic scaling of the distance is the 
best one can get for homological codes based on a closed 
(two-dimensional) surfaces. For such codes he proved that 

< c(log fc)^n, (1) 

with a constant c. Hyperbolic surface codes could thus offer 
the possibility of encoding more logical qubits for a given 
number of physical qubits and a desired level of protection or 
distance. 

However, these codes arise from tilings of a hyperbolic 
surface, endowed with a hyperbolic distance function, lead¬ 
ing to a two-dimensional lay-out of qubits in which qubits 
interacting in parity check measurements are 0(1) distance 
away with respect to this hyperbolic distance function. If 
one embeds such a tiling of the surface in an ordinary two- 
dimensional plane this embedding cannot preserve the distance 
function (i.e. cannot be isometric). This means that some qubits 
interacting in parity check measurements will be long-distance 
as compared to interactions between other qubits. This non¬ 
isometric embedding is demonstrated in the Poincare disk 
model of the hyperbolic plane shown for a specific tiling in 
Fig.ffl One way to visualize the hyperbolic plane (see more 
details in Section IIII-Cl l is by mapping all of its points into 
the interior of a unit disc (radius R = V) within the euclidean 
plane. The shortest line through two given points, a geodesic, is 
mapped to an arc of a circle on which both points lie and which 
intersects the boundary of the disc orthogonally. This mapping 
preserves angles (it is conformal) but it heavily distorts areas 
and distances. For example, all pentagons in Fig. [T] have 
the same hyperbolic area while clearly in the Poincare disk 
model the distances become exponentially smaller around the 
boundary (Other models of the hyperbolic plane exist such 
as the Poincare half-space model in which one represents 
= {x,y G M I y > 0}, i.e. a half-plane, with a metric 
ds^ = {dx^ + dy‘^)/y‘^). 

In the quantum code constructions below, qubits are placed 
at the edges of the tiles. The tiles themselves represent parity 
Z-checks while tiles on the dual lattice (star operators centered 
at a vertex) represent X-parity checks of the code. In order to 
encode (multiple) qubits, the hyperbolic surface has to be 
closed, similar as a torus is obtained by closing the euclidean 
plane by taking periodic boundaries. For the hyperbolic plane 
H^, if we want this closing procedure to be isometric, so that 
locally the surface looks like H^, it implies that one obtains a 
torus with many handles (or genus g > 1), with each additional 
handle encoding 2 additional qubits (see the GauB-Bonnet 
theorem in Section IIII-Db . Of course, one can also make a 
many-handled torus out of a euclidean plane, but this would 
require some stretching or crumpling; this is another way of 
saying that the flat curvature of the plane cannot be maintained 
in the many-handled torus. 

For practical applications it is more interesting to have a 
fully planar 2D qubit layout such as in the surface code than 


a code which is defined on qubits distributed over the two- 
dimensional surface of a many-handled torus. For this reason 
we consider how to planarize the hyperbolic surface codes 
in Section IIV-CI we will show that the creation of explicit 
boundaries on which logical operators have to start and end 
reduces the efficiency of these codes so that a planar code with 
constant rate k/n> ci necessarily has a constant distance, no 
longer increasing (logarithmically) with n. 

In order to construct a hyperbolic surface code one can 
pursue the following strategy. First, one constructs a tiling of 
the hyperbolic plane: such a tiling should be isometric in that 
it respects this distance measure of the hyperbolic plane, see 
Section im] Such tiling can be associated with a group such 
that group elements label elementary triangular tiles in which 
the basic tiles can be decomposed. In order to compactify the 
hyperbolic plane, i.e. make it into a multi-handled torus which 
encodes some qubits, one identifies elementary triangular tiles 
which differ by certain translations: this can be done by 
finding a normal subgroup of the tiling group. Finding all 
possible normal subgroups (which normal subgroup one picks 
determines how many qubits k are encoded) is not simple 
and has to rely on previous mathematical constructions. In 
m, im regular hyperbolic tilings were considered and a 
construction of finding certain normal subgroups due to Siran 
was used. In a regular tiling {r, s} each tile is a regular r- 
gon and s regular r-gons meet at each vertex in the tiling (see 
Section Imll. Ref. an] considered the {7,3} hyperbolic tiling 
and studied numerically the distance of the logical operators 
for different encoding rates k/n. 

Our work expands on these previous constructions by (1) 
using mathematical tools that allow one to find any compacti- 
fication of any (r, sj hyperbolic tiling, (2) analyzing the noise 
threshold of several interesting tilings and (3) considering the 
logical error probability and distance of encoded qubits of 
interesting tilings. Asymptotic scaling of [[n,/c,(i]] codes is 
interesting, but in a physical implementation one is interested 
in using a code with some fixed parameters. We show how 
a hyperbolic surface code can outperform the regular surface 
code in terms of encoding more qubits with the same number 
of physical qubits for a fixed distance (see Table HHl in Section 
HV^ . 

For superconducting transmon qubits, the interaction dis¬ 
tance between transmon qubits is flexible when the interqubit 
coupling is mediated via a transmission line (as in the ar¬ 
chitecture described in ifThl or the experiments in Q, ll8l). 
The interaction length-scale can vary between 100 ym (for a 
compact bus-resonator) to 10 cm for a transmission line, hence 
a factor 10^. It may be of interest to consider how to convert 
the graph generated by the hyperbolic tiling into a planar grid- 
graph or a double-layer grid-graph with equidistant qubits and 
minimal non-planar, crossing edges. 

In a recent paper Ell the authors introduce holographic 
quantum codes on which physical qubits live at the endpoints 
of a hyperbolically-tiled surface (for the holographic pentagon 
code defined in that paper, this is the {5,4} tiling) such 
that these degrees of freedom represent the encoding of 
logical degrees of freedom which are explicitly entered 
through the bulk. These codes are thus not directly related 
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to the hyperbolic surface codes analyzed in this paper in 
which the physical qubits occupy the entire hyperbolic surface. 

We will review the definition and decoding of homolog¬ 
ical codes in Section El In Section [HI] we discuss regular 
tilings and the construction of regular tiled, closed surfaces. 
We first consider the euclidean case and then generalize to 
manifolds which admit constant curvature. In Section ng 
we discuss concrete examples of quantum codes based on 
closed orientable hyperbolic surfaces with regular tilings. In 
Section IIV-BI we give a numerical analysis of the thresholds 
of families of hyperbolic surface codes. Finally, we discuss 
the construction of hyperbolic surface codes which are planar 
(they are embeddable on a disc). 

II. Homological codes 

In this section we review how a tiled, closed surface can 
be used to dehne a homological (CSS) quantum code using 
Z 2 -homology. 

A. 'Ij 2 -homology 

In Z 2 -homology, statements about the topology of a surface 
are turned into statements of linear algebra over the held Z 2 = 
{0,1} in which all operations are carried out modulo 2. 

Consider a tiled, closed surface X. We assume that two faces 
overlap on at most a single edge and two edges overlap on at 
most a single vertex. From now on we will call faces, edges 
and vertices 2-cells, 1-cells and 0-cells. The tiled surface X is 
called a cell complex. 

For a given cell complex X the subsets of i-cells in X form 
a Z 2 -vector space where the addition of two subsets is given 
by their symmetric difference (the set of elements which are 
contained in one set or the other but not in both). The standard 
basis of this space is given by all sets containing a single i- 
cell. We will identify these sets with the standard basis of 
Z™‘ where rrii is the number of i-cells contained in X. We 
denote these vector spaces as Ci{X) or simply Ci and call 
their elements i-chains. For i = 1, 2 one can dehne a boundary 
operator di 

di '. Ci ^ Ci—i ( 2 ) 

by its action on the basis vectors; let e® G be a single z-cell, 
then di{C) G Ci-i is the sum of all {i — l)-cells incident to 

e®. 

An important property of the boundary operator is that a 
boundary does not have a boundary 

81082 = 0, (3) 

which is equivalent to the statement that boundaries of faces 
are closed loops. One can also dehne the dual map called the 
coboundary operator 

(4) 


which assigns to each i-cell all (j-l-l)-cells which are incident 
to itOl Viewed as matrices, the coboundary operator is the 
transpose of the boundary operator of one dimension higher 

= (5) 

One uses both mappings to dehne subspaces of Ci. From the 
boundary operator one obtains the cycle space Z\ = ker 81 and 
the boundary space Bi = im 92- Due to Eqn. we always 
have the boundary spaces as subspaces of the cycle spaces 
Bi < Zi. 

Intuitively, the cycle space Z\ represents all 1-chains which 
have no boundary, i.e. collections of closed loops. Certainly, 
all boundaries of 2-chains have no boundary (which is the 
content of Eqn. O but there might be loops which “go around” 
the surface in a homologically non-trivial way. We can count 
the latter ones by introducing the first homology group Hi — 
ZijBi, such that dimiFi is the number of loops that are non¬ 
trivial. Eor an orientable surface, dimiJi is twice the number 
of its “handles” (genus) g of the surface. 

The coboundary operator defines the cocycle space Z^ = 
kerdi and the coboundary space B^ = im Sq with B^ < Z^. 
Erom Eqn. |5] it follows that 

6)ci+i = (a, 9i+i6)ci (6) 

where {■,-)ci denotes the standard inner product in Ci. 
Together with Eqn. |3] it follows that the cocycle space is 
the orthogonal complement of the boundary space and the 
coboundary space is the orthogonal complement of the cycle 
space, i.e. 

B^ = Z^ and Z^ = B^t- (7) 

In other words, the coboundaries are exactly those chains 
which have even overlap with the cycles and the cocycles 
are exactly those chains which have even overlap with the 
boundaries. 

The i-cells of a tiled surface X correspond to (2 — i)-cells 
of the dual tiling X*. By linear extension, this gives rise to an 
isomorphism 

* : CfiX) ^ C2-.(X*) (8) 

of the i-chains of X to the (2 — i)-chains of X*. We hrst note 
that going to the dual chains leaves the inner product (even or 
oddness of their overlap) invariant 

{a,b)c, = {*a,*b)c‘_.- (9) 

Directly from their definitions, it is also clear that applying the 
coboundary operator to a chain of a cell complex is equivalent 
to going to the dual complex and applying the boundary 
operator of the complementary dimension 

Si = o 82 -i o * (10) 


*In this paper we identify chains with cochains via the inner product. 


Si : Ci ^ Ci+i 
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Or in diagrammatic form; 


C, -^ Q+i 

* * 

^2-1 <-^2-i-l 


( 11 ) 


Eqn. [To] allows us to interpret the coboundaries and 
cocycles as boundaries and cycles of the dual structure 
(Fig.ia. 



Fig. 2: Part of the {6,3}-tiling. The dual tiling (dashed) is the 
{3,6}-tiling. The boundary of a face (red) corresponds to the 
coboundary of a vertex in the dual structure (blue). 


B. Quantum codes 

From the Z 2 -homology of a tiled surface, we can dehne a so 
called stabilizer code encoding k logical qubits into n physical 
qubits, TL = (C^)®" (see background on stabilizer codes in 
01). The 2^-dimensional codespace C Q B is defined as the 
+1 eigenspace of an Abelian subgroup S {—I ^ S) of the 
Pauli group. To turn a closed, tiled surface X into a stabilizer 
code we identify all its edges with qubits. The boundaries of 
the faces are used to define Z-type check operators (stabilizer 
elements). For every face we add a generator to S which 
acts as Z on all edges which belong to the boundary of 
the face. Equivalently, the set of coboundaries of vertices 
gives a generating set of 2f-check operators. Note that in 
general the generating set of Z- and AT-stabilizers is not 
independent (Fig. O. From now on, when we refer to the 
stabilizer generators we mean this canonical set associated to 
the faces and vertices. 

We call the number of physical qubits on which a Pauli 
operator acts non-trivially its weight. If the surface is tiled 
by a regular tiling given by the Schlafli symbol {r, s} (see 
definition in Section ITlI-Ab . then r is the weight of the Z- and 
s the weight of the 2f-stabilizers. 

The number of physical qubits is simply the number of edges 
n = dim Ci. The number of encoded qubits can be calculated 


Fig. 3; A {4,4}-tiling of the torus (opposite sides are identi¬ 
fied). The 1-chain highlighted in black is the boundary of the 
gray square, but also the boundary of all white squares. 


by taking the number of physical qubits and subtracting the 
number of restrictions that the stabilizers impose: 

k = dim Cl — dimi?^ — dimi?i 
= dim Cl — dim Z^ — dim Bi 
= dim Cl — (dim Ci — dim Zi ) — dim Bi 
= dimTTi 

There is another way of seeing this which helps the intuitive 
understanding of homological codes; Consider the subgroup 
N{S) of the Pauli group which consists of those operators 
which commute with all elements in 5'0 The action of N{S) 
on T-L leaves the space of logical qubits as a whole invariant. 
The Z- and AT-type elements of N{S) stand in one-to-one 
correspondence with the cycles and cocycles by Eqn. |7| All 
stabilizers (the boundaries and coboundaries) have a trivial 
action on the logical qubits. The operators which are cycles 
but not boundaries, i.e. non-trivial loops, have a non-trivial 
action on the qubits. With the help of Eqn. |5] one can see 
that cycles which are not boundaries and cocycles which are 
not coboundaries come in dimi/i pairs with odd overlap, 
whereas the overlap of elements from different pairs is even. 
We identify these pairs of non-trivial loops as the X- and Z- 
operators, each acting on one of the logical qubits. 

The distance d of the code is given by the minimum weight 
of a logical operator. By the previous discussion this is the 
same as the minimum length of a homologically non-trivial 
cycle in the tiling or its dual. This quantity is also known as 
the combinatorial systole and denoted by csys (or csys* for 
the dual case). 


C. Error correction in homological codes 

A stabilizer code is used to actively correct errors by 
measuring the generators of S. The measurement result is 
called the syndrome. The syndrome information is then used 
to infer a recovery operation. Note that the error only needs 
to be corrected up to a stabilizer element. 

For homological codes, error correction can be understood 
geometrically. Every Pauli error can be identihed with two 
chains Ex,Ez G Ci, where the support of each chain tells 

^The N stands for normalize!'. As Pauli operators either commute or anti¬ 
commute N{S) coincides with all elements that commute with S. 
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US where the error acts as X respectively Z on the qubits. 
Consider for concreteness Bz- The syndrome is the set of all 
X-stabilizers which anti-commute with the error Bz (Fig. |5]l. 
Since we identify X-stabilizers with vertices the syndrome is 
simply the support of dBz- Analogously, the syndrome of the 
X-error Bx is given by SBx or equivalently by the boundary 
of *Bx in the dual. 

A Z-error Bz is corrected by applying a Z-type Pauli 
operator with support R G Ci such that its boundary coincides 
with the syndrome. After this is done the system is back in a 
code state. If the sum of chains Bz + i? is in i?i, we applied 
a stabilizer and no operation was performed on the logical 
qubits. If, however, Bz + R contains a non-trivial loop, then 
we applied a logical operator and our encoded information is 
corrupted. An X-error is corrected in the exact same way but 
in the dual tiling. 

In general, a syndrome will consist of many vertices of 
the (dual) lattice. The problem of correcting the error is 
choosing and connecting pairs such that the likelihood of a 
resulting logical operator is minimized. This depends on the 
probabilities for different errors to occur. Those are in turn 
determined by the error model. Typically, errors with high 
weight are assigned a low probability. Hence, a good recovery 
strategy is to find the lowest-weight chain R consistent with 
the syndrome. This problem is a variant of minimum weight 
(perfect) matching (MWM) which can be solved efficiently by 
the so-called blossom algorithm. 

III. Tilings of closed surfaces 
A. Euclidean tilings 

In this section we introduce tilings of the euclidean plane 
(denoted by E^). The concepts introduced in this section will 
directly carry over to spaces with constant curvature. 

The set of all distance-preserving mappings from the eu¬ 
clidean plane onto itself forms a group under composition. This 
group is called the group of isometries Isom(E^). Its elements 
are either translations, rotations or combinations of translations 
and reflections (called glide reflections). Let G be a subgroup 
of Isom(E^). For a point P GEf, we call the set of all points 
that P is mapped to by elements of G the orbit of P under 
G. 


Stabilizer code 

Homology 

Sz 

Bi 

Sx 

or Bl 

N{Sx)z 

Zi 

N(Sz)x 

Z^ or Z* 

n 

dimGi 

k 

dim Bti 

d 

min(csys, csys*) 


Fig. 4; Overview: Corresponding notions in the language of 
stabilizer codes and homology. The subscript X and Z restricts 
the set to X- respectively Z-type operators. 



Fig. 5: The error chain Bz is highlighted. At its boundary we 
have two vertices which constitute the syndrome of the error. 


Take for example the group of all integer translations to 
the right. The orbit of a point P = {x,y) G £2 is the set 
{{x + n,y) I n G Z}. 

A fundamental domain of a group G is a part of the plane 
which contains a representative of each G-orbit with at most 
one representative of each G-orbit in its interior. For the 
aforementioned example a choice for a fundamental domain 
is the strip [0,1] x R. 

The action of G on the points of induces an action on 
a fundamental domain P of G. We denote the application of 
a g G G on P hy B * g. The identity element e G G acts 
as the identity map P * e = P and the action is compatible 
with the group multiplication. By the latter we mean that the 
application of gh on P is the same as first applying g to P 
and then h to the result, or in short P * (gh) = {P * g) * h. 

If we have a set of (finite) parts of the plane which cover 
E^ and all parts either do not overlap or overlap on their 
boundaries we call this a tiling of E^. Each element of the 
tiling is called a tile or face. The orbit of a fundamental 
domain naturally gives rise to a tiling. 

We will focus on regular tilings where all faces are (con¬ 
gruent) regular polygons which are placed edge-to-edge. Each 
regular tiling can be labeled by the number of sides of the 
polygons r and the number of polygons s meeting at a point. 
This label is known as the Schldfli symbol {r, s}. Every {r, s}- 
tiling has a well defined dual tiling with Schlafli symbol {s, r} 
(Fig.©. 

We will now discuss Wythoff’s kaleidoscopic construction 
which allows us to relate a regular tiling {r, s} to a group of 
isometries G = Gr,s- 

A regular r-gon has 2r symmetries generated by the re¬ 
flections on its symmetry axes (Fig. j^. The symmetry axes 
induce a triangulation of the r-gon into 2r (right) triangles. 
The triangulation of the faces induces a triangulation of the 
whole tiling into triangles with internal angles 7r/2,7r/r and 
tt/s. 

Let Gr,s be generated by the reflections on the edges a, b, c 
of a single triangle. We assume that a, b and c are arranged in 
clockwise order such that a is opposite to the angle •n/r and 
b is opposite to the angle tt/s. 

A reflection applied twice is the same as doing nothing. 
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SO that a? = iP' = (? = e. Additionally, two consequtive 
reflections on two lines that intersect at an angle 7r/fc for 
some positive integer k correspond to a rotation around the 
intersection point by an angle 27r/fc. Thus, k rotations give 
the identity. 

These are all the relations that a, b and c fulfill and the 
multiplication of any two elements of Gr,s is completely 
determined by them. Hence, we can write compactly 

Gr s = (a, b, c I = b^ = = (ab)^ = (be)’" = (ca)’’ = e). 

( 13 ) 

This is called the presentation of the group Gr,s- One way 
to think about group presentations is as a set of strings or 
words consisting of {a,b,c,a~^,b~^,c~^}. Group multipli¬ 
cation simply corresponds to a concatenation of words. All 
words which differ by a subword that is equivalent to e by the 
rules of group multiplication (such as g~^g — e) or any of 
the given relations are considered to be equal. Abstract group 
presentations can be directly used in computer algebra systems 
such as GAP or Magma. 

There is an important subgroup of Gr,s which we will 
denote by G+j,. It consists of all orientation-preserving maps 
(rigid motions). The group G+g is generated by all double 
reflections p = be and a = ca, so that 

Gt,s = (P, cr I p’' = cr" = (per)^ = e). (14) 

The generators p and a act as a clockwise 27r/r and 27r/s 
rotations, respectively. Note that the G+^-orbit of a triangle 
only covers “half” the plane as we cannot map between 
triangles that share an edge. We fix this by taking two triangles 
related by a 6-reflection as the fundamental domain of G+g. 

The important observation to make is that (by construction) 
each element of the orbit of the group acting on a fundamental 
domain is uniquely labeled by a group element. By fixing a 
triangle Fq that is the fundamental domain of G+g, every other 
triangle F can be labeled by the group element g G G+j, that 
maps Fq onto it, i.e. Fo*g = F. This allows us to forget about 
the group action (Fq was arbitrary in the first place) and only 
talk about G+g. 

Let (p) = {e, p, p^,..., p’'“^} be the cyclic subgroup of 
G+s that is generated by the rotation p. The faces of the {r, s}- 
tiling (the r-gons) are in one-to-one correspondence to the 
triangles up to a rotation by p. In group theoretic language 
these are the left cosets of the subgroup (p) denoted by g{p) = 
{PjPP, pp^,... ,gp’'~^} for a p G G+^. Similarly the vertices 
and edges of the tiling can be uniquely labeled by left cosets 
of the cyclic subgroups (tr) and (per), respectively. 

Note that faces, edges and vertices are incident if and only 
if their associated cosets share a common element. This means 
that the topology of the lattice is encoded in the group Gf^ 
as well. 

B. Quotient surfaces 

In this section we review how to use isometries to define 
closed, two-dimensional manifolds (or surfaces) which locally 
look like E^. 



(a) Action of the symmetry group on a single face. 



(b) Rotations acting on the lattice 


Fig. 6: Group acting on the {6,3}-tiling. 


Let H he a subgroup of isometries acting on E^. We can 
construct a new space E^/iJ where we identify all points that 
differ by the application of an element of H. Formally, each 
point is now a set of the form {P * g \ g G H}, i.e. the orbit 
of some point P. 

To avoid degenerate cases, one demands that the action of 
F[ does not have any fixed points. Additionally, we assume 
that there exists an e > 0 such that for all points P in E^ the 
discs of radius e around every point in the orbit P * H do not 
overlap. If H meets these requirements one says that it acts 
fixed-point free and discontinuous. In this case there exists a 
fundamental domain of H around every point of E^. 

In an earlier example we considered a group F[ consisting 
of all integer translations along the x-axis. H acts fixed-point 
free and discontinuous and its quotient surface E^/iT is an 
infinite cylinder with circumference one. One can imagine that 
H wraps E^ up in its x-direction and positions all points 
belonging to one orbit ‘on top of each other’. In this sense 
E^ is a covering of the surface E^/iJ and H is called the 
covering group. 

Note that in the example the geometry of the cylinder is 
the same as the geometry of the plane within discs of radius 
< 1/2 . This radius is called the injectivity radius Rinj. The 
injectivity radius plays an important role in code construction 
of Section [H] as it provides a lower-bound on the length of the 
shortest closed loop on 'Ef /H that cannot be contracted to a 
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point. 

For a quantum code it follows that if the tiling is defined on 
a quotient surface then all loops contained within a disc with 
radius Rmj must be boundaries. This is because the geometry 
within the disc is the same as the geometry of and all loops 
in are boundaries of the set of faces that they surround. The 
injectivity radius therefore provides a lower bound on the code 
distance d, assuming that all edges are of unit length. 

Which quotient surfaces E^/iF admit a regular tiling? One 
condition on H is, that it needs to respect the tiling structure. 
This means that it must be a subgroup of the tiling group Gr,s- 

The faces, edges and vertices of the tiling on the quotient 
surface E^/iF are similarly labeled by the action of Gr,s on 
an arbitrarily chosen face. However, for this to work on the 
quotient surface it must “look the same” everywhere and in 
every direction. The action of the covering group H should 
not affect the labels. This is the case if for all h G 7T and 
g G Gr,s we have ghg~^ G H. We say that iF is a normal 
subgroup of Gr,s- 

If H does not contain any glide-reflections (H is a subgroup 
of G+s), its quotient surface ¥?/H is orientable. 

In the tiling of E^/iJ, each face can be labeled by a set 
of the form g{p)H = {gp^h \ n G {1,... ,r}, h G H} for 
a p G G^g. Similarly, we have a labeling of the edges and 
vertices of the quotient surface E^/iF using cosets of {pa)H 
and {a)H. Faces, edges and vertices of E^/iJ are incident if 
and only if their associated cosets share a common element. 

The number of cosets of a subgroup in a parent group is 
called the index of a subgroup. If the index of H in Gr,s is 
finite then E^/iJ has finite area. 

In conclusion, one is able to define a group of isometries 
Gr,s which encodes a regular tiling {r, s}. The problem of 
finding closed, euclidean, orientable surfaces which admit a 
{r, s} tiling reduces to finding finite index, normal subgroups 
of Gr,s which have no fixed-points. 

C. Surfaces with Curvature 

The curvature of an arc of radius R can be measured by 
1/R. We can likewise measure the curvature of a surface at 
a point P by cutting it with two normal planes to create two 
arcs with radii Ri and i? 2 - For a sphere, the two resulting arcs 
have their radii of curvature on the same side. For a saddle, the 
radii of curvature can show in different directions. In this case 
we assign the radii different signs. We define the curvature 
K at point P by taking the maximum and minimum possible 
values of i?i and R 2 and setting K = Hence, the 

curvature of a sphere with radius Ris K = 1/R^ at any point. 
For a saddle, K will be negative and for the euclidean plane 
K vanishes. K is called the Gaussian curvature. 

Any triangle on the unit sphere has angles which sum up 
to more than tt due to the curvature. This is called an angular 
defect. We can define a regular tiling on using Wythoff’s 
construction that we introduced in Section IIII-AI A regular r- 
gon on can be divided into 2r triangles. For a regular tiling 
{r, s} the internal angles of the triangles are 7r/2, 7 r/r and 
tt/s. To make these angles to add up to more than tt we have 
the condition that l/r-|-l/s > 1/2. The only regular tilings 


which fullflll this condition are {3,3}, {3,4}, {4,3}, {3,5} 
and {5,3}. These tilings are well known as they correspond 
to the Platonic solids. 

Similarly, for the euclidean plane we have the condition that 
for a {r, s}-tiling it has to hold that 1/r -f 1 /s = 1/2 as all 
angles need to sum up to tt. Hence, {6,3}, {3,6} and {4,4} 
are the only regular tilings of E^. 

If the curvature is negative the internal angles have to sum 
up to a value smaller than tt. The condition l/r-fl/s < 1/2 
is now met for an infinite number of pairs r, s. The plane with 
constant, negative curvature K = —1 is called the hyperbolic 
plane (see Fig. 7). 



Fig. 7: The {7,3}-tiled hyperbolic plane in the Poincare 
disc model. The fundamental domains of Gt^s are filled in 
black and white. Fundamental domains of the same color are 
related by an element of G 7 3 . 

All constructions introduced in Section IIII-AI can also be 
applied to the sphere and the hyperbolic plane H^. We 
then obtain (regularly tiled) quotient surfaces with curvature 
K — +1 and K = —1, respectively. 

Curvature does not only determine what regular tilings are 
possible but it also determines the topology of the quotient 
surfaces. The sphere itself only has contractible loops. The 
only quotient surface it admits is called the projective plane, 
which is non-orientable. It is obtained by factoring out a 
reflection across a plane which intersects the sphere at a great 
circle. Homological quantum codes derived from tilings of the 
projective plane were analyzed in ||4|. 

Since the euclidean plane has no curvature, all translations 
commute. Therefore two non-colinear translations suffice to 
make the resulting quotient surface closed. If the covering 
group H contains translations only, the surface V?/PI is a 
torus. If iJ also contains a glide-reflection, the quotient surface 
E^/iJ is a Klein bottle. These are the only two closed surfaces 
with vanishing curvature. 
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We will see in this section that the hyperbolic plane admits 
more freedom and allows for quotient spaces which have high 
genus. 

We end this section by mentioning that (up to scaling) all 
closed surfaces of constant curvature can be expressed as 
quotient surfaces V?/H or 'M?/H according to the 

Killing-Hopf theorem M- 


D. Properties of hyperbolic surfaces codes 

In this section we will analyze the general properties of 
closed, hyperbolic surfaces and the quantum codes obtained 
from them. 

The Gaufi-Bonnet theorem connects the curvature of sur¬ 
faces to their topology. It states that for a closed, oriented 
surface M the integral over the Gaussian curvature K is 
proportional to the genus g of M 


[ KdA = 2TT{2-2g). 

JM 


(15) 


Hence, negative curvature of M, K = — 1 at any point enforces 
a proportionality between the genus and the area; 


2 g = 


area(M) 

27r 


+ 2. 


(16) 


It is well-known im that the area of a regular, hyperbolic 
r-gon n is given by its angular defect 


area(n) = (r — 2 ) 7 r — ^ Qfi (17) 

i=l 


where ai are the internal angles of H. If M is tiled by {r, s} 
we have ai = 27r/s and thus 


area(M) = ttF - 2 - 2 ^) (18) 

where F is the number of faces in the tiling of M. 

By substituting Eqn. [18] in Eqn. [T 6 | we see that a [[n, k, dj] 
quantum code derived from the {r, s}-tiling of M will have a 
rate 



(19) 


The asymptotic rate (n oo) only depends on the tiling and 
the rate is higher for larger values of r and s. In the previous 
section we established that there are infinitely many hyperbolic 
tilings subject to the constraint l/r -f 1/s < 1/2. The tiling 
giving the smallest combined weight r -f s of the stabilizer 
generators is the {5,4}-tiling (and, of course, its dual {4, 5}). 
Compare the above result to the rate of codes obtained from 
tilings of oriented quotient surfaces of the euclidean plane. As 
discussed in the previous section they are topologically all tori. 
Hence, they only encode two logical qubits, regardless of the 
number of physical qubits. 

What about the distance of hyperbolic surface codes? Unfor¬ 
tunately, there exists no formula for the combinatorial systole 
of a tiled quotient surface. 

It has been proven in l(20ll that for a {r, s}-tiling of 
for any given g G N there exists a hyperbolic surface H^/id 


tiled by {r, s} with combinatorial systole larger than q and 
n < C"^s/2, where n is the number of edges and 

( 20 ) 

This means that there exists a family of quantum codes 
where the distance is lower bounded by a function growing 
logarithmically in the number of physical qubits (see also 
mi Since the logarithm is base C and the value of C is very 
larg 43 , this bound is not relevant for any practical purposes. 

Unfortunately, we cannot do better than a logarithmic 
growth of the distance. In ll22ll the author proves that for a 
{r, s}-tiled quotient surface the combinatorial systole d has an 
upper bound 

d<^^og^{ 2 n). ( 21 ) 

One can also compare this to the bound in Eqn. [T] 

We can mention that hyperbolic surface codes of type 
{3,oo}, where infinitely many triangles meet at so-called ide¬ 
alized vertices at infinity, have been studied in ll2^ . However, 
these surfaces have a dual systole csys* G 0(1), hence a 
distance upper bounded by a constant Il24l . 

Once a hyperbolic code is constructed, we can efficiently 
compute its distance using the following algorithm due to 
Bravyi 123: Eirst we compute a basis of the cocycle space 
= ker(5i) via GauBian elimination on the coboundary 
operator (5i, which takes 0 {'nf ) steps where n is the number of 
edges. Remember that the cocycles correspond to the logical 
X operators (see Eig. |4|. 

Second, we consider two copies of the lattice X and X'0 
Eor every element h in the basis of we define a new lattice 
D using the following procedure: If the edge e = {u, v) in 
X and e' = (u', v') in X' are in the support of b they are 
removed from X and X'. At the same time we introduce new 
edges connecting u in X with v' in X' and vice versa u' 
with V. We basically take the original lattice X and paste in 
a copy X' along the coboundary b to obtain D. A non-trivial 
cycle has to have odd overlap with a cocycle, as established 
in Section III-BI Hence, the length of the shortest non-trivial 
cycle in X which intersects b and passes through point v can be 
determined by computing the shortest path between v and v' in 
D. This can be done using, for example, Dijkstra’s algorithm 
which takes 0{n log(n)) steps. By looping over the k elements 
of the cocycle basis b and base points v we can thus determine 
the systole using 0{n^ -|- kn'^ log(n)) steps. The cosystole is 
determined by repeating this procedure using the dual lattice. 

IV. Constructions 
A. Examples of Hyperbolic Surface Codes 

In Section IIII-BI we established a correspondence between 
{r, s}-tiled, oriented quotient surfaces and fixed-point free, 
finite index, normal subgroups of the group Gf^. In ll26l the 

^There exist a follow-up result 1211 which involves number theoretic 
functions and gives a C of some magnitudes smaller. 

^Since we are only interested in paths on the lattice we can consider them 
to be graphs and forget about the faces in either lattice. 
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n 

k 

d 

csys 

csys’' 

translation 

48 

4 

3 

6 

3 


168 

14 

4 

8 

4 


384 

32 

4 

12 

4 

Gp-^T 

648 

54 

6 

14 

6 

(Tp‘*‘ap ^(jp^a ^ p^ (J ^p ^<yp ^ 

768 

64 

6 

16 

6 

Z t Z — i\>5 B — 1 — B — Z 

Gp \p G ) p G p Gp 


TABLE II; Examples of compactified {8, 3}-lattices. The com¬ 
binatorial systole of the primal lattices is larger than that of 
the dual lattices. 


authors present an algorithm that, given the presentation of 
G+j, (Eqn. fT4l i. finds all of its normal subgroups iT up to a 
given index. 

In Table |I] and Table El we present some examples of 
possible compactifications. The first three columns contain 
the three code parameters, followed by the combinatorial 
systole of the primal and the dual lattices. The last column 
contains translations which compactify the lattice. In the group 
theoretic language of Section IIII-BI the translations and all 
conjugations by elements of Gr,s generate the group H. If the 
compactification is achieved by a single translation, then the 
length of this translation determines the combinatorial systole. 

B. Numerical Analysis of Codes 

In this section we present the results of simulated error 
cotTections on hyperbolic surface codes. Our error model is 
that of independent X and Z errors in which a qubit can 
undergo independently an X error with probability p and 
a Z error with probability p at each time-step. After these 
errors happen one applies the MWM-decoder which tries to 
infer what error occurred. Hereby we assume that syndrome 
measurements can be done perfectly. The decoder succeeds if 
the product of the real and the inferred error is in the image of 
the boundary operator. We do this independently for the dual 
lattice as well. 

We gather statistics by repeating this procedure N times. 
The probability of a logical error Piog on any of the encoded 
qubits is estimated by taking the ratio between failed correc¬ 
tions and number of trials. In the usual simulation of the toric 
code one plots the logical error rate against one type of error 
{X or Z) only, as the quantum error correction for both types 
of errors is identical. This is not the case for general {r, s} 
codes. 

The results of a simulation with N = 4x 10"* trials are shown 
in Eigure[8] Eor increasing probability of a physical error p the 
probability for a logical error occuring on any of the qubits 
Piog approaches 1 — (1/2)^*^. This is because the decoder can 
do no worse than random guessing and the probability to guess 
correctly for all 2k logical operators Aj, Zi is (1/2)^^ (for the 
toric code this is 0.75 for a single type of Pauli error). 

In a numerical simulation of the noise threshold of the toric 
code (see e.g. 1221) the logical error probability curves for 
different numbers of qubits n all intersect at one point, namely 
the threshold pc, clearly pointing to the fact that the noise 
threshold corresponds to a phase-transition. In Eig. we see 
that the distributions shift to the left with increasing number 



physical error rate p 


(a) The zoomed-in plot shows p in the interval [0.005, 0.025] and 
Piog in the interval [0,0.06]. The black line in the zoomed-in 
plot represents p = Piog- 



physical error rate p 


(b) Piog-p-diagmm for {5,4}-tiled surfaces. The examples were 
chosen such that the ratio of distance to physical qubits forms a 
strictly increasing sequence. 


Eig. 8; Logical error Piog against physical error p for {5,4}- 
tiled surfaces. Every data point was obtained from 4 x 10'^ 
runs. 
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n 

k 

d 

csys 

csys’' 

translation(s) 

60 

8 

4 

6 

4 

{{ap-^Yp-^r 

160 

18 

6 

8 

6 

ap {ap ) p <y p <yp 

360 

38 

8 

8 

8 

(crp-‘<jY{p-^a-^p-Y‘‘ 

1800 

182 

10 

10 

10 


1920 

194 

10 

12 

10 



TABLE I; Examples of compactified {5,4}-lattices. The examples were chosen such that we have a sequence with growing n 
and d. The second to last lattice is compactified by two independent translations. 


of qubits. Note that even if the distance is not increasing we 
can still observe a better protection against errors. We believe 
that this is the case due to an increase in the average minimum 
weight of the operators of the logical qubits. 

When taking a sequence of {5,4}-codes with strictly in¬ 
creasing distance and number of physical qubits, we see 
that the first three lines cross (Eig. ISFt . However, the line 
corresponding to the {5,4}-code with n = 1800 qubits does 
not cross. The data is somewhat inconclusive in showing 
whether there is a noise threshold and what its asymptotic 
value is. Eurther examples can be found in the Appendix lAl 

It is important to note that it has been shown in ll28ll that 
there exists a finite threshold for all homological codes with 
distance lower bounded by a logarithm in the number of qubits, 
when quantum error correction is perfect or noisy. 

Eor hyperbolic surface codes, we expect the existence of 
another threshold where all logical qubits are potentially cor¬ 
rupted. In f2^ the authors analyze the two distinct perculation 
thresholds for a prototype model of a hyperbolic lattice. 
Here we only consider the threshold where an arbitrary qubit 
becomes corrupted. 

Our data suggests that for the {5,4}-tiling there exists a 
threshold between 1% and 5% (see Eigure EJl. Eor the toric 
code the threshold against this error model is pc = 10.3% as 
determined in Ez). 

The distance of the code depends on the parameters r and 
s of the tiling. Eor increasing r, s (and thus a code with better 
rate) the lower bound on the distance becomes smaller and one 
would expect that the threshold also goes down (at least the 
lower bound on the threshold in EHl becomes smaller). 

In our simulations of various {r, s} codes we see that the 
threshold does go down for increasing stabilizer weight r 
and s which is conform to our expectation (see Eigure [T3] in 
Appendix IATi. 

To compare the toric code to the hyperbolic surface codes 
we consider the overhead to produce logical qubits and protect 
them against decoherence. To do this we mark up to which 
physical error rate a code can protect all of its qubits with 
probability at least 0.999. The results are shown in Eigure 0 

We see that the toric codes protect against errors up to a 
higher physical error rate. However, their overhead is quite 
big compared to the hyperbolic surface codes. 

Eor each family of codes the encoding rate decreases with n 
(Ean.fTgll. Since we are below threshold and Piog goes to zero 
for increasing n the capability of the codes to protect against 
errors increases. Eor each code family we obtain a slope which 
converges to the asymptotic rate of the code. Codes with lower 


bD 

C 

■v 



physical error rate p 


Eig. 9: Encoding rate of a code which protects qubits with 
probability > 0.999. The number of physical qubits varies 
between 60 and 960. Data points are labeled by the tiling. The 
two instances of the toric code are L = 4 and L = 6. 


rate offer better protection against errors. Since the weight 
of the logicals is not a monotone function in the number of 
physical qubits the slopes are not monotone either. 

The hyperbolic surface code with the smallest encoding rate 
is the {7, 3}-code which was analyzed in ifTSl . There are only 
two {7, 3} codes with less than 10^ qubits. The next smallest 
rate is achieved by the {5,4}-code which appears to offer the 
best error protection of all hyperbolic surface codes. 

In Table [HI] we show the minimum weight of the logical 
operators of some hyperbolic surface codes. The weights of 
the logicals were determined by brute force as the method to 
determine the distance from Section IIII-DI does not directly 
apply. 

The smallest instance of the {5,4}-code encodes k = 8 
logical qubits into n = 60 physical qubits and has distance 
d = A. An efficient use of the surface code (in which the lattice 
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Lattice 

n 

k 

d 

weight of Z-logicals 

weight of X-logicals 

{5,4} 

60 

8 

4 

6(8) 

4(8) 

{5.5} 

30 

8 

3 

3 (6), 4 (2) 

3 (8) 


40 

10 

4 

4(10) 

4 (9), 5, (1) 

{6,4} 

36 

8 

4 

4 (4), 6 (4) 

4(8) 

{6,5} 

60 

18 

3 

4 (15), 6 (3) 

3, (3), 4, (13), 5 (1), 6 (1) 

{6.6} 

54 

20 

4 

4 (18), 6 (2) 

4(20) 


60 

22 

4 

4 (15), 5 (6), 6 (1) 

4 (16), 5 (6) 


TABLE III: Minimum weight of the logical operators of some 
small hyperbolic surface codes. The number in parenthesis 
is the number of distinct logical operators with the same 
minimum weight. 


is chopped off at the boundaries, see lIMIl ) has parameters [[n = 
cP,k = l,d]] (instead of [[n = cf + {d — = l,(i]])- 

Hence to encode 8 qubits with the surface code with distance 
4 requires 128 physical qubits, showing that the hyperbolic 
construction can lead to more efficient coding than the surface 
code. 

Note that the surface code has a higher threshold, making 
it more favourable for noisy qubits. For more coherent qubits 
hyperbolic surface codes might offer an advantage. 

C. Hyperbolic surface codes with open boundaries 

One can try to construct a hyperbolic surface code with 
open boundaries in various ways. Such codes, if they have 
good distance and rate, would have a larger practical appeal 
than the hyperbolic surface codes, as they would correspond 
to planar graphs. 

A hrst idea is to cut open the multi-handled torus to obtain a 
surface with some boundaries. One method produces a surface 
with punctures, the other method would produce a region as 
depicted in Fig. (TO] 

In m it was first described how a surface with punctures 
can encode multiple qubits. This type of encoding is used in 
the surface code architecture described in IISTlI where qubits 
are encoded in both rough or smooth double-holes punctured 
in the surface. For this encoding which uses both smooth and 
rough holes, there is a simple scheme with which one can 
perform a CNOT gate, first described in ll32]| . 

One could image obtaining such a punctured surface by 
cutting open the multi-handled torus. One would start with a 
torus with g handles and no boundary, nor punctures, encoding 
k = 2g qubits. If one cuts open one handle, one has two 
punctures (and g — 1 handles) in the remaining surface. Cutting 
open another disjoint handle similarly generates two punctures 
and cutting open the last handle leads to a surface with 2g—l 
punctures and one outer boundary (the last puncture). For a 
tiled surface one can make a smooth cm or a rough cut. For a 
smooth cut, one follows a non-trivial Z-loop 7 on the graph 
around the handle and all edges e G on this loop 7 are 
replaced by double pairs of edges, i.e. Ej —> E:^ U Ef^ so 
that the two Z-plaquettes adjoining these edges either act on 
edges in Ei^ or E!^ (but not both). At each vertex on this path 
the number of A-check operators is doubled such that one X- 
check operator acts on the edges in E^ and the other on the 
edges in E^. Note that in this procedure, one adds as many 


qubits as one adds A-check operators as the number of edges 
and vertices around a loop are the same. For the many-handled 
torus, there is one linear dependency between all plaquette 
Z-checks (the product of all them is /) and similarly one 
linear dependency for the A-checks. For the final punctured 
surface (with smooth holes), there is only a linear dependency 
between all A-checks: this directly implies that the punctured 
surface will encode 2^ — 1 logical qubits. Another way of 
seeing that a surface with a smooth outer boundary and 2g—l 
smooth punctures or holes encodes 2^ — 1 logical qubits is 
by enumerating the logical operators: for each hole the logical 
Z is a Z-loop around the hole, while the logical A is a A- 
loop on the dual lattice to the outer boundary (see e.g. m). 
However, one does not need to use all these logical qubits. 

In the smooth double-hole encoding described in ED, ID, 
one uses a pair of holes to encode one qubit where the logical 
Z is the Z-\oop around any of the two holes and the logical 
A is the A-distance between the two holes. In this way, the 
surface with 2^—1 holes can encode g double-hole qubits 
where the last double-hole qubit is formed by the last hole 
plus the boundary ‘hole’. 

However, does this procedure preserve the distance of the 
original code and is it even possible to execute this procedure 
on the graph obtained from a {r, s}-tiling of the closed 
hyperbolic surface? The answer is in fact no. If we have a 
family of {r, s} codes with increasing k ^ n and distance 
lower bounded as some c(r, s) log n, then logical Z loops 
(which act on at least on c(r,s) log n qubits) must overlap 
on many qubits. If there are H(/c) such non-intersecting loops, 
then there must be at least f2(fclogn) = r2(nlogn) qubits 
in total, which is not the case! Hence, one cannot find such 
a set of non-overlapping loops along which to cut as logical 
operators must share a lot of space (support). 

Another method of cutting the many-handled tiled torus to 
create smooth and rough boundaries would allow one to create 
a coding region as in Fig. [TO] One can encode multiple qubits 
into such a surface code El by using alternating rough and 
smooth boundary of sufficiently large length. The boundary of 
the encoding region is divid^ into 2k regions and encodes k — 
1 logical qubits. The logical Z operators can start and terminate 
at rough boundary regions while the logical A operators start 
and terminate at smooth boundary regions. One can consider 
the asymptotic scaling of the parameters n,k and d of the 
code represented by such an encoding region, where we do 
not assume anything about the distribution of qubits inside the 
region, meaning that for such a class of codes with increasing n 
qubits can get closer together (or further apart) if we represent 
the qubits on the euclidean plane. One can simply argue that 
for such a family of homological surface codes with boundaries 
the following bound should hold 

kd < cn, (22) 

with a constant c. This bound shows that for codes based on 
tiling a surface with open boundaries, it is not possible to 
have a constant encoding rate k/n and a distance increasing 
as logn. We argue this bound for the encoding region in 
Fig. [TOl but similar arguments should hold for the tiling of 
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a surface with holes. The simple argument goes as follows. In 
order to encode k — 1 qubits, one divides the set of boundary 
edges inbound into 2k regions. Clearly, the total number of 
qubits n > inbound- The logical operators run from rough to a 
rough or smooth to smooth boundary, hence their distance is 
upper-bounded by the length of a region c|inbound |/fc < cn/k 
for some constant c. This is true as one can always let the 
logical operator run along (or close to) the boundary edges. 
This results in Eqn. |22] If the encoding region is a disk with 
many holes, we can argue similarly. We enumerate the number 
of edges (qubits) around the non-trivial disjoint holes which 
should be less than the total number of qubits n. The number 
of logical qubits k scales linearly with the number of holes 
and the distance d is at most the number of edges around each 
hole, resulting in Eq. |22] 

Thus we may ask whether using hyperbolic geometry gives 
any advantage over euclidean geometry when considering 
tiled surfaces with a boundary. The previously proven bound 
kd^ < cn is still worse than the simple bound kd < cn. It is 
worth noting that encoding via the boundaries as in Eig. [10] 
has worse scaling than an optimal surface code encoding in 
the case of a euclidean metric on the underlying qubits. In 
this case, the perimeter of the polygon scales at most as y/n 
where n is the total number of qubits. This implies that the 
distance of the code d < Cy/n/k, which is worse than the 
kd? < cn bound. If we use hyperbolic geometry and imagine 
the encoding region as a partially-tiled hyperbolic plane, see 
the explicit construction below, then the number of qubits at the 
boundary scales like the total number of qubits, which seems 
promising. However, the minimum weight logical operators 
which run from boundary to boundary, will run along shortest 
paths, geodesics, which go through the interior of the region. 
We show a small explicit example of such a code in the next 
section, but the underlying construction is unlikely to give good 
asymptotic scaling behavior. 

D. Small Planar Code Example 

We generate a planar graph with a hnite number of vertices 
by taking a regular tiling {r, s}, partially tiling the hyperbolic 
plane and then modifying the boundary so that the planar 
graph encodes logical qubits, as in Eig. do] An example of 
the procedure to create such a graph encoding multiple qubits 
consists of the following steps; we illustrate the idea in Eig.fTTI 
for the {5, 5} tiling; 

1. Eor a {r, s} tiling, start with a single r-gon (call it level-1 
r-gon) and reflect in the edges of this r-gon to generate 
level-2 r-gons. Repeat this for the level-2 r-gons etc. so that 
one obtains a planar graph G = {V,E) where the faces are 
the level-1 to level-m r-gons for some m. In the Eigure, we 
have started with 4 level-1 r-gons and generated only the 
level-2 r-gons and have stopped there. With every face / we 
associate a plaquette Z-check acting on the qubits on the 
boundary of this face, and with every vertex v of this graph 
one associates a star -check acting on the qubits on edges 
adjacent to the vertex. At the boundaries the weight of the 
AT-checks will be two or s (for even s) or s — 1 (for odd s), in 
the interior the X-checks have weight s. The code associated 



Eig. 10; A polygon with 2k sides, alternatingly ‘rough’ and 
‘smooth’ (see definitions in Section lIV-Dl i. can encode k — 1 
logiral qubits a. Shown is fc = 4 and a choice for the logical 
Xi, Zi operators. 


with this starting graph G encodes no qubits. There is one 
linear dependency between all X-checks and so the number 
of linearly independent AT-checks is V — 1. There is no linear 
dependency between the Z-checks of which there are F. As 
this is a planar graph (with Euler characteristic x = 1)> ons 
has E = V — 1 + F where E is the number of edges, equal 
to the number of qubits n, hence no encoded qubits. 

All boundaries in this graph are so-called smooth boundaries 
at which a string of X-errors can start and end. Eor this graph 
which encodes no qubits, such a string can be annihilated 
by star X-checks, hence there are no logical operators. Thus 
we need to modify this graph in order to create so-called 
rough boundaries. At rough boundaries a string of Z-errors 
can start. If one alternates rough and smooth boundaries as 
in Eig. [T0| these strings can no longer be annihilated by X- 
or Z-checks, but have to run from boundary to boundary. 
What is the procedure for creating several rough boundaries 
at which Z-strings can end? 

2. Eirst, for the given starting graph G, one counts the 
number of boundary edges in E'bound as | inbound |: these are 
defined to be the edges that the level-fc r-gons would be 
reflected in to generate level k + 1 r-gons. In order to create a 
code where all encoded qubits have about the same distance, 
one wants to divide the set E'bound in 2fc equally-sized sets or 
regions with fc > 2 so that each subset of size [|£'bound|/(2/c)J 
corresponds to a rough or smooth boimdary. Such code can 
encode k — 1 qubits IS); the logical Zi, i = 1,... ,k — 1 
will run from the ith to the (i + l)th rough boundary, while 
the logical Xi will run from the smooth boundary in between 
the 2 th and {i -f l)th rough boundary to the {i + l)th smooth 
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boundary (see Fig. fTOl l. If k is chosen too large, then each 
smooth/rough region would be too short and the shortest path 
between two regions would not go through the bulk. Hence one 
should choose k such that the shortest path along the perimeter 
is about the same length as the shortest path through the bulk; 
in that case one uses the available qubit-space optimally. This 
choice is illustrated in Fig. [TT] where one has 60 edges along 
the boundary which we divide up into 10 = 2/c regions of 6 
edges each, thus encoding 4 logical qubits. 

The creation of a rough boundary in a region of edges 
^bound — -E^bound consists of the following 3 steps (variants 
are possible). First, one removes all X-check operators of 
weight-2 which have 2 edges among F^^ou'nd weight-2 

X-check has an edge in one also removes it). One 

could in principle also remove checks with weight more than 
2 at the boundary but this makes the lattice a bit smaller, so 
in this construction we prefer to remove only the weight-2 
checks. Then all qubits on which only a single plaquette 
Z-check acts are removed from the lattice and these plaquettes 
are thus modified (in the toric code they are the weight-3 
plaquettes at the boundary). The removal of the X-checks 
makes it possible for a Z-string to start on an edge attached to 
a removed X-check. We need to make sure that such a string 
can only run from rough boundary to another rough boundary. 
Since we did not remove some X-checks in the rough region 
(some weight-5 checks in FigiriTl. in general of weight s), 
one needs to add weight-2 ZZ checks to the stabilizer, as 
these operators commute with all the current checks. These 
weight-2 checks are indicated (in red) in Fig. [TT] For the 
construction in^ig. [TT| one can verify that the minimum 
weight logical Z is of weight-4, while the minimum-weight 
logical X is of weight-5; we draw examples of these logical 
operators in the Figure. The total number of physical qubits is 
n = 65. The most efficient use of the surface code (in which 
the lattice is chopped off at the boundaries, see EOl) has 
parameters [[cP, 1, d]] (instead of [[d^ + {d— 1)^, 1, d]]). Hence 
to encode 4 qubits with the surface code with distance 4 
requires 64 physical qubits and distance 5 would require 100 
qubits, showing that this simple hyperbolic construction can 
lead to a somewhat more efficient coding than the surface code. 

In order to construct a planar graph encoding k — 1 logical 
qubits, one can also generate the planar starting graph G 
encoding no qubits by repeatedly rotating an elementary 
r-gon by 27r/s around its vertices. By this rotation, one 
generates a new generation of r-gons, the level-2 r-gons etc. 
In this construction, one removes all X-checks in the region 
where the rough boundary has to be formed. This means that 
no weight-2 ZZ checks need to be added. Only vertices on 
which a single plaquette Z-operator acts are to be further 
removed. 


V. Conclusion 

We have discussed explicit constructions of families of 
2D topological quantum codes derived from regular tilings 
on surfaces with negative curvature. On a hyperbolic surface 



Fig. 11; A [[65,4,4]] code based on the {5,5} tiling. The 
distance_of the logical X is in fact 5 while it is 4 for the 
logical Z. The number of boundary edges in inbound of the 
starting graph was 60 and was divided into 10 regions each 
with 6 edges. Shown in red are the additional weight-2 ZZ 
checks. 


without boundaries these hyperbolic surface codes have an 
asymptotically constant rate and low 0(1) weight stabilizer 
checks and distance scaling as logn for a total of n physical 
qubits. This trade-off of distance versus encoding rate is much 
more favorable than for the toric or surface code. For practical 
purposes it is desirable to minimize the weight of the check 
operators as the quantum circuit through which one measures 
the check operators is noisy. For this reason, the most attractive 
code may be the {5,4} hyperbolic surface code with an 
asymptotic rate of k/n ^ 1/10; we have obtained some partial 
evidence that the noise threshold of this code for noiseless eiTor 
correction is above 2%. It may of be interest to consider the 
explicit usage of these codes for a 2D quantum memory given 
that they allow for a higher-qubit storage density. 
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Appendix 



physical error rate p 
(a) {8, 3}-tiling 



physical error rate p 
(b) {6, 6}-tiling 



physical error rate p 
(a) {7, 3}-tiling 



physical error rate p 
(b) {6, 4}-tiling 


Fig. 13: Logical error Piog against physical error p for hyper¬ 
bolic surface codes with higher weight. We chose the examples 
to be strictly increasing in the distance. Each data point was 
obtained by W = 4 x 10^ trials. 


Fig. 12: Logical error Piog against physical error p for hyper¬ 
bolic surface codes with higher weight. Each data point was 
obtained by W = 4 x 10^ trials. 


































